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Abstract 

For a recurrent linear diffusion on R+ we study the asymptotics of 
the distribution of its local time at as the time parameter tends to 
infinity. Under the assumption that the Levy measure of the inverse 
local time is subexponential this distribution behaves asymtotically as 
a multiple of the Levy measure. Using spectral representations we find 
the exact value of the multiple. For this we also need a result on the 
asymptotic behavior of the convolution of a subexponential distribu- 
tion and an arbitrary distribution on R+. The exact knowledge of the 
asymptotic behavior of the distribution of the local time allows us to 
analyze the process derived via a penalization procedure with the local 
time. This result generalizes the penalizations obtained in Roynette, 
Vallois and Yor [22] for Bessel processes. 
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1 Introduction 



1. Let X be a linear regular recurrent diffusion taking values in R_|_ with 
an instantaneously reflecting boundary and +oo a natural boundary. Let 
Px and Ej; denote, respectively, the probability measure and the expectation 
associated with X when started from x > 0. We assume that X is defined in 
the canonical space C of continuous functions u : R+ i— > R+. Let 

Ct := a{u;{s) : s < t} 

denote the smallest cr-algebra making the co-ordinate mappings up to time t 
measurable and take C to be the smallest cr-algebra including all cr-algebras 
Ct, t > 0. 

We let m and S denote the speed measure and the scale function of X, 
respectively. We normalize S by S{0) — and remark that 5'(+cxd) = +00 
since we assume X to be recurrent. It is also assumed that m does not have 
atoms. Recall that X has a jointly continuous transition density p{t;x,y) 
with respect to m, i.e., 



,{Xt e A)^ J^p{t;x,y)m{dy), 



where A is a Borel subset of R_|_. Moreover, p is symmetric in x and y, that 
is, p{t; X, y) = pit] y, x). The Green or the resolvent kernel of X is defined for 
A > via 

poo 

Rxix,y) / e-^'p{t;x,y)dt, (1.1) 
Let {Lf^^ : ^ > 0} denote the local time of X at y normalized via 

LJ^) = lim ^ , / l[y,y+5){^s) ds. (1.2) 



sio m{{y,y + S)) Jq 

For y = we write simply Lt, and define for £ > 

n := inf{s : > £}, (1.3) 

i.e., T :— {t£ : i > 0} is the right continuous inverse of {L^}. As is well known 
T is an increasing Levy process, in other words, a subordinator and its Levy 
exponent is given by 

Eo (exp(-AT^)) = exp {-e/RxiO, 0)) 

POO 

= exp(-£ / u{dv){l - e-^"')), (1-4) 
Jo 
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where z/ is the Levy measure of r. The assumption that the speed measure 
does not have an atom at implies that r does not have a drift. 



2. We are interested in the asymptotic behavior of the distribution of Lt as 
t tends to infinity. The basic assumption under which this study is done 
is the subexponentiality of the Levy measure of r (see Section HI). The 
subexponentiality assumption is equivalent with the relation (cf. Proposition 

P{n>t) ~ £z/((t,+oo)) V£>0. 
Here and throughout the paper the notation 

f{x) ~ g{x), 

x—i-a 

where / and g are real valued functions and a is allowed to take also "values" 
+00 or —00, means that 

i„„44=i, 

x^a g[x) 

Since r is the inverse of L, it also holds (see Proposition 14. 10 

Po(^*<^) ~ ^v{{t,+^)). 

t— >+oo 

To extend this for an arbitrary starting state x > 0, we first show that (see 
Proposition 14. 2p 

V^{H,>t) ~ 5(x)z/((t,+oo)), 
where Hq := inf{t : Xt = y}, and then (see Proposition 14. 3h 

PxiLt<£) ~ iS{x)+i)iy{{t,+oo)). (1.5) 

>+oo 

Our motivation for relation (11.51) arose from the desire to generalize the 
penalization result obtained for Bessel processes in Roynette, Vallois and Yor 
[22J (see also [19j and l2lj). From our point of view, since many of the pe- 
nalization results are derived for Brownian motion and Bessel processes, it is 
important to increase understanding of the assumptions needed to guaran- 
tee the validity of such results for more general diffusions. In particular, we 
prove that (see Theorem 15.21 and Example 15. Sp 

MHLt) \C ) ^ g^xMLu) + 1 - HiL^) =: a.s., (1.6) 
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where h is a probability density function on R+ (with some nice properties) 
and H is the corresponding distribution function. 

3. The paper is organised as follows. In the next section basic properties 
on subexponentiality are presented and a new result (Lemma 12. 4p on the 
limiting behavior of the convolution of an subexponential and a more general 
distribution is derived. In Section 3 we study the spectral representations 
of the hitting time distributions and the Levy measure. In Section 4 results 
on subexponentiality and the spectral representations are combined to yield 
relation (jl.Sp . Hereby we also need a weak form of a Tauberian theorem given 
as Lemma EH] in Appendix. The application in penalizations is discussed in 
Section 5. To make the paper more readable we state and prove first the 
general theorem on penalizations. After this the penalization with local time 
is treated and (II. 6p is proved. The paper is concluded by characterizing the 
law of the canonical process under the penalized measure induced by the 
martingale M^. Using absolute continuity and the compensation formula for 
excursions we are able to shorten the proof when compared with the one in 



2 Subexponentiality 

In this section we present some basic results on subexponential probability 
distributions. Later, in Section [H it is assumed that the probability dis- 
tribution induced by the tail of the Levy measure of r is subexponential. 
This assumption allows us to deduce the crucial limiting behavior of the first 
hitting time distribution (see Proposition 14. 2p . 

Definition 2.1. The probability distribution function F on (0, +cxd) such 
that 

F(0+) = 0, F{x) < 1 Va; > 0, lim F{x) = 1 (2.1) 

x—*oo 

is called subexponential if 

lim FVF{x) /F{x) = 2 (2.2) 

X— > + oo 

where * denotes the convolution and F{x) := 1 — F{x) the complementary 
distribution function. 
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For the following two lemmas and their proofs we refer Chistyakov [3] and 
Embrechts et al. [5]. 

Lemma 2.2. If F is a probability distribution function satisfying 112. 1\) and 

F{x) ~ x-'^H{x) 

with a > and H a slowly varying function then F is subexponential. 
Lemma 2.3. If F is subexponential then 
(i) uniformly on compact y-sets 

lim + = 1, (2.3) 

(a) for all £ > 0, 

lim e^^F(a;) = +oo (2.4) 

The proof of the next lemma uses some ideas from Teugels [25] p. 1006. 

Lemma 2.4. Let F and G be two probability distributions on R+. Assume 
that 

(1) F is subexponential, 

(2) lim,^ooG(x)/F(a;) = c>0. 
Then 

lim YVGix)/ (G(x) + F{x)) = 1. (2.5) 

Proof. Let e G (0, 1). By assumption (2) there exists 6 = 6{e) such that for 

X > 6 

c (1 - e)F{x) < G{x) < c (1 + e)F{x). (2.6) 

Observe that 



F * G{x) = 1- F* G{x) = 1 - F{x) + F{x) - G{x - y) dF{y) 

Jo 

= F{x)+ [ G{x-y)dF{y). 
Jo 
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We assume now, throughout the proof, that x > S and write 



G{x) + F{x) G{x) + F{x) G{x) + F{x) 

where 

G(x-y) , 
h{x) := / dF{y) 



and 



Gix) 



l^-s G{x) 



Obviously, by assumption (2), the claim (12.5p follows if we show that 

lim Ii{x) = 1 

X— >oo 

and 



lim l2{x) = 0. 
Proof of i2.9\) . Since G{x — < 1 we have 

rfF(y) F{x)-F{x-S) 



h x) < 



G{x) G{x) 

F{x -6)- F{x) 



G{x) 
F{x) fF{x - 6) 



1 . 



G{x) V F{x) 

Using now (12.3^ and assumption (2) yields (I2.9p . 
Proof of (ESj. Since _ _ 

G(x - > Gix) 

we have 

Jo 6- a; 



Consequently, 



liminf Ii{x) > 1. 

z— >+oo 
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To derive an upper estimate, notice first that 

h{x) < —- / -^=7^ dF{y), (2.11) 
l-e Jo F{x) 

because, from ( I2.6p . 

x>5 G{x) >c{l- e)F{x) 

and 

y < X — 5 =^ x — y>6 =^ G{x — y) < c{l + e)F{x — y). 
Next we develop the integral term in fl2.1ip as follows 

-<5_ 

Fix - y) dF{y) 

-s 

{l-F{x-y)) dF{y) 

F{x-6)- / F{x-y)dF{y) 
Jo 

px—5 

F{x) - F{x- y) dF{y) + F{x -5)- F{x) 
Jo 

F{x)~ r F{x-y)dF{y)- r dF{y) 

Jo Jx-6 

px~5 

F{x)- / F{x-y)dF{y) 
Jo 

X PX 

F{x - y) dF{y) - F{x - y) dF{y) 

-S Jx~S 

X 



PX 

< F{x) - / F{x -y)dF{y). 
Jo 



Hence, 



x—5 



F{x - y) dF{y) < F{x) -F* F{x) = F* F{x) - F{x). 
Consequently, from ( I2.1ip . 

hix) < 



l+e\ f F*F{x) - F{x) 



1-eJ \ F{x) 
7 



and using (I2.2p and letting e ^ we obtain 

limsup Ii{x) < 1 

x— >+oo 

which together with (12.10p proves (I2.8P completing the proof of Lemma 12.41 

□ 



3 Spectral representations 

Spectral representations play a crucial role in our study of asympotic proper- 
ties of the hitting time distributions. In this section we recall basic properties 
of these representations and derive some useful estimates. For references on 
spectral theory of strings, we list [7], [9], [4], [TO], fT4], fT2] , p], fU], and 

m- 

Besides the diffusion X itself, it is important to study X when killed at 
the first hitting time of 0, denoted X = {Xt '■ t > 0}, i.e., the diffusion with 
the sample paths 

X,:=|^*' (3.1) 
\d, t>H„ ^ ' 

where Hq := inf {t : Xf = y}, and (9 is a point isolated from R+ (a "cemetary" 
point). Then {Xt : t > 0} is a diffusion with the same scale and speed as X. 
Let p denote the transition density of X with respect to m : 

PxiXt G dy) = P^iXt edy-t< Hq) = Pit; x, y) m{dy). (3.2) 

Recall that the density of the Pa;-distribution of Hq exists and is given by 

U{t) := P.iHo e dt)/dt = hm^^l^. (3.3) 

Moreover, the Levy measure z/ of the inverse local time r, see (II. 2p and 
(11.31) . is absolutely continuous with respect to the Lebesgue measure, and 
the density of z/ satisfies 

.(AO/*. = Um^ (3.4) 
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We define now the basic eigenfunctions A{x;'y) and C{x;'j) associated 
with X and X, respectively, via the integral equations (recall that S is con- 
tinuous and m has no atoms) 

A{x--f) = l-j dS{y) / m{dz)A{z;-f), 
Jo Jo 

C{x-^)=S{x)-^ r dS{y) r m{dz)C{z-^), (3.5) 
Jo Jo 

and the initial values 

^(0;7) = 1, A'(0;7):=lim^^^^7\^ = 0, (3.6) 

xiO 

C(0;7) = 0, C'(0;7):=nm^^ = l. (3.7) 

xio b [x) 

Let {Anjand {C„} be two families of functions defined by 



Ao(x) = l, An+i{x)= dS{y) / m{dz)Aniz) (3.8) 

Jo Jo 

and 

Co{x) = S{x), C„+i(x)= rdS{y) r m{dz)C^{z), (3.9) 

Jo Jo 

respectively. Then the functions A(a;; 7) and C(x;7) are explicitly given by 

00 

A{x-^) = Y,{-irM^)- (3.10) 

n=0 

and 

00 

C{x-^) = Y,{'irCn{x), (3.11) 

n=0 

respectively (see Kac and Krein [7J p. 29). In the next lemma we give an 
estimate which shows that the series for C converges rapidly for all values 
on 7 and x > 0. A similar estimate for A can be found in Dym and McKean 
[4] p. 162. 
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Lemma 3.1. The functions x ^ C„(x), x > 0, n = 0, 1, 2, . . . , are positive, 
increasing and satisfy 

Cn{x) < ^ Six) M{y) dSiy)^ (3.12) 

where M{z) = m(0, z). 

Proof. The fact that C„ are positive and increasing is immediate from (13.9P . 
Clearly (13.121) holds for n = 0. Hence, consider 

Cn+i{x) = / dS{y) / m{du)Cn{u) 
Jo Jo 

< f dS{y) r m{du) — S{u)( f M{z)dS{z) 
Jo Jo ^- \Jo 

<^S{x) [ dS{y) r m{du) ( [ M{z)dS{z) 
Jo Jo \Jo 

<-^(a;) dS{y) M{z)dS{z)] M{y) 

n+l 

S{x) ( / M{y)dS{y) ' 



(n + l)! 

where we have used the facts that x ^-^ S{x) is increasing and x i— > M{x) is 
positive. □ 

Lemma 3.2. The function x ^— ^ C{x; 7) satisfies the inequality 

|C(a;;7)| < S{x) exp ( \-f\ [ M{z)dS{z)] . (3.13) 



Proof. This follows readily from (13.111) and ( I3.12p . □ 

From Krein's theory of strings it is known (see |4] p. 176, and [7l [T4l [T3]) 
that there exists a a-finite measure denoted A, called the principal spectral 
measure of X, with the property 

A(dz) 

-^—^ < 00 3.14 

z+1 ^ ' 
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such that the transition density of X can be represented as 

POO 

Pit; x,y)= / e-^* 7) Aiy; 7) A(t/7). (3.15) 



We remark that from the assumption that m does not have an atom at it 
follows (see [4j p.l92) that A([0, 00)) = 00. 

Analogously, for the killed process X there exists (see [T2], [E]) a cr-finite 
measure, denoted A and called the principal spectral measure of X, such 
that 

and ^ 

A(dz) 

-^-^ = 00. (3.17) 
The transition density of X can be represented as 







00 



Pit; x,y)= e-^* C(x; 7) Ciy; 7) A(rf7). (3.18) 



The result of the next proposition can be found also in [15]. Since the 
proof in [1^ is not complete in all details we found it worthwhile to give here 
a new proof. 

Proposition 3.3. (i) The density of the P^- distribution of the first hitting 
time Hq has the spectral representation 

POO 

/.o(t)= / e^^*C(x;7)A(c?7). (3.19) 



(a) The density of the Levy measure of the inverse local time at has the 
spectral representation 

POO 

u(t) = / e-^*A(rf7). (3.20) 
Jo 

Proof, (i) Combining f l3.3p and (13. 18^ yields 

r /,x 1. P{t;x,y) 
fxoit) = hm— -— . 
yio S{y) 

= lim re-^*C(x;7)^^A(rf7). 
yio Jo S{y} 
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We show that the limit can be taken inside the integral by the Lebesgue 
dominated convergence theorem. Let t > be fixed an choose e such that 



t - / M{z)dS{z) > t/2. 



Then, from Lemma [3.21 for 7 > and < y < e we have 

^^^^ < exp (^-7 - j'^ M{z)dS{z)j^ < e'^^/^ 
Consequently, it remains to show that 

/■oo 

/ e-^*/2 \C{x] 7) I A(rf7) < 00. (3.21) 
By the Cauchy-Schwartz inequality 

e-^*/2|C(x;7)|A(rf7)l 

/•OO noo 

< / e-^*/2(C(x;7))'A(d7) / e-^'/^A{d^) 
Jo Jo 

poo 

= p{t/2;x,x) / e-^*/2A(rf7). 
Jo 

Clearly, p{t/2;x,x) < 00 and, by flXTB]) . e-^^/^A{d-f) < 00. These es- 
timates allow us to use the Lebesgue dominated convergence theorem and 
since (cf. i KTh ) 

limC(y;7)/5(y) = C'(0;7) = l 

y^O 

the proof of (i) is complete. Representation (13.20p can be proved similarly 
using formula (13.41) . ( 13.19^ . (13. 7p and the estimates derived above. We leave 
the details to the reader. □ 

Remark 3.4. Consider 

/•oo /"OO 

[lAt)u{t)dt= / dt{lM) A(c/7)e-^* 



A(ci7) / dt{lM)e-^\ 
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A straightforward integration yields 

/•OO 1 

/ (lAt)e-^*(it = — (l-e-T), 
Jo 7 

and, consequently, {\3.16i) is equivalent with (cf. [T^) 

/"OO 

/ (1 A t) z>(t) cit < OO, 
Jo 

which is the crucial property of the Levy measure of a subordinator. For 
(El^j, see p. 82. and fl^ . 

Example 3.5. Let R = {Rt : t > 0} and R = {Rt : t > 0} be Bessel 
processes of dimension < 5 < 2 reflected at and killed at 0, respectively. 
We compute explicit spectral representations associated with R and R. 

From, e.g., p. 133 the following information concerning R and R can 
be found: 

Speed measure 

m(dx) = 2x^-2° cix a:={2-5)/2. (3.22) 

Scale function 

5(x) = ;^x2" (3.23) 
Transition density of R (w.r.t. m) 

p{t- X, y) = ^^{xyr exp (^-^^±^^ (^^) , x,y > 0. (3.24) 

Transition density of R (w.r.t. m) 

p{t; X, y) = ^^{xyr exp la (y ) , x,y > 0. (3.25) 

To find the Krein measure A associated with R we exploit formulas 113.1^) 
and Ili3. 24\ ) with x = y = and use 
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to obtain 

^—{l-a) poo 



p{t; 0, 0) = lim p(t; x, y) = ——— ^ = / e"^* A{d^). 







Inverting the Laplace transform yields 

We apply formula 1^3. 1(h) . 113. 22) . and ^3. 23\) to find the function A{x; , and, 
hence, compute first directly via ^3. 8\) 

AJx) = — ^ , / -, n = 0,1,2,.... 

^ ' 2"r(n + l)r(n+l-a) 

Consequently, after some manipulations, we have 

A{x; 7) = r(l - a) 2-- (xv^)" (^v^) , 
where J denotes the usual Bess el function of the first kind, i.e., 

w ^ _ V (-l)"(^/2)-+^- 

Z.r(n + l)r(z. + n+l)' 

and, finally, putting pieces together into 113.1^) yields 

p{t;x,y) = ^ ^ e-^'{xyrj^^ (a^v^) J-a (vV^) d^. (3.27) 

Next we compute the Krein measure A associated with R. For this, we deduce 
from (E3j, (Q, ( TOP , and ^3M) 

p(t;x,y) 2i-"ar(^+") /""^ -:^ . , x 
z> t = lim ^ \ = -— =/ e-^*A rf7 , 3.28 

and, consequently, inverting the Laplace transform gives 
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Similarly as above, we apply formula ^3.11\) to find the function C {x; 'j) , and, 
hence, compute first directly via ^3. 9\) 

CJx) = , ^ \ , -, n = 0,1,2,.... 

Consequently, after some manipulations, 

C{x; 7) = T{a) 2("-')/2 (^^^^^ . 

and 

Pit; x,y) = ^ e-^* {xy^ Ja (a; v^) J. (y v^) ^7- (3.30) 
5*66 also Karlin and Taylor p. 338. 

Example 3.6. Taking above a = 1/2 yields formulas for Brownian motion. 
Recall 

[Y [Y 

Jiniz) = \ — sinz, and J-i/2{z) = \ — cos 2;. 
Consequently, from ^3.27\) 

POO 1 

p{t; x,y) = - e-^* cos(a; v^) cos{y^) (3.31) 
Jo V ^7 

^ ^g-(x~y)V(2i) ^ g-(x+y)2/(2t) 



2V27rt 
and from ^3. 3(A) 

1 f°° 
P{t;x,y) = - e 
Jo 



7i sm{xy/2^) sm{y^) ^ji^^^ 



2^j2^t 

Moreover, 



27 



1 I -{x-yfl{2t) _ (x+y)V(2t) 



1 /"°° T 

/.o(t) = - / e-^* sin(xv/2^) = ——= e'^^^^^ 

TT Jo t^/^\/2TT 



Ht) = - / e-^* v^ci7 = 7^777^- (3.32) 
71^ Jo t-^/^v27r 



and 



From ( foil) we obtain A{d-f) = d-f/{^^^/2^), and from ( TOl) A(d7) = 
v/27rf7/7r. 5*66 ako ifar/in and Taylor ^ p. 557 and 393, and p. 120. 
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Proposition 3.7. (i) The complementary P x-distribution function of Hq has 
the spectral representation 

Px{Ho>t)= / -e-^*C(x;7)A(d7). (3.33) 
Jo 7 

(a) The Levy measure has the spectral representation 

poo /»oo -| 

z/((t,oo))=/ u{s)ds= -e-^*A(rf7). (3.34) 
Jt Jo 7 

Proof. Formulas (13.33^ and ( I3.34p follow from ( 13.19^ and ( I3.20p , respectively, 
using Fubini's theorem. To obtain (13.34^ is straightforward but for (13.330 the 
applicability of Fubini's theorem needs to be justified. Indeed, from (13.19^ 
we have informally 

PxiHo >t)= fxois) ds= ds A(c?7) e-^^ C(x; 7) 
Jt Jt Jo 

POO 

A(d7) / dse-^'C{x;-f) 



leading to (13.331) . To make this rigorous, we verify that for all a; > 

f 00 -1 

-e-^*|C(x;7)| A(d7) < 00. 
7 

Consider first for e > 



Ki:= r-e-^*|C(x;7)|A(rf7). 
Jo 7 



lo 7 

By the basic estimate ( I3.13p for < 7 < e 



|C(x;7)| < S(x) exp J M{z)dS{z) 
and, consequently, 

Ki < S{x) exp (^e M{z)dS{z)^ ^ e"^* A(rf7) < 00 
by (I3T611 . Next, let 

i^2:= / -e-^*|C(^;7)|A(d7). 

J e 7 
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By the Cauchy-Schwartz inequality 



/oo n 



oo 

e-^*(C(x;7))'A(rf7). 



The first term on the right hand side is finite by (13. 161) . For the second term 
we have 



e-^*(C(x;7))'A(rf7)< / e"^* 7))' A(rf7). 

< p{t; X, x) < 00. 



The proof of f l3.33p is now complete. □ 

4 Asymtotic behavior of the distribution of Lf 

as t ^ +(X) 

We make the following assumption concerning the Levy measure of the in- 
verse local time process {ti : i > 0} valid throughout the rest of the paper 
(if nothing else is stated) 

(A) The probability distribution function 

xf-^— — — -, X > 1, 
1/(1, +00) 

is assumed to be subexponential. 

It is known, see Sato [24] p. 164, that Assumption (A) is equivalent with 

P{re>t) ~ £z/((t,+cx))) V^>0, (4.1) 

and also with 

The law of Tg is subexponential for every i > 0. (4.2) 
Proposition 4.1. For any fixed i > 0, it holds 

Po{Lt<i) ~ iiy{{t,+^)). (4.3) 

t—f + OO 
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Proof. The claim follows immediately from (14.11) since 

Po{Lt <i) = P(r, > t). 

□ 

Our goal is to study the asymptotic behavior of Lt under P^;. For this, we 
analyze first the distribution of the hitting time Hq. The proof of the next 
proposition is based on Lemma [6T] stated and proved in Section [6] below. 

Proposition 4.2. For any x > 0, it holds 

P,{Ho>t) ~ 5(x)z/((t,+oo)). (4.4) 

Proof. Recall from (13.33^ and (13.34^ in Proposition 13.71 the spectral repre- 
sentations 

/"°° 1 

P.,{Ho >t)= - e-^* C{x- 7) A(rf7) (4.5) 
Jo 7 

and 

zy((t,+oo))= / -e-^*A(rf7). (4.6) 
Jo 7 

We apply Lemma EH with /i((i7) = A((i7)/7, (71(7) = C(x;7) and (72(7) = 
S{x). Then, the mapping t Px{Hq > t) has the role of fi and t \—y 
S{x) z/((t, +cxd)) the role of /2. Condition (16.11) takes the form 

lim S{x)u{{t,+oo))e^' = 

t^oo 

and this holds by Assumption (A) and (I2.4p . Moreover, condition (16. 2p means 
now 

lim C(a;; 7)/S'(a;) = 1 
and this is true since using estimate (13.121) in (13. lip we obtain 

< a|7|e'^l^l 



C{x;^) ^ 



Six) 

with some a and (3 depending only on x. Consequently, (16.3P in Lemma [6. II 
holds and, hence, the proof of the proposition is complete. □ 

The main result of this section is as follows. 
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Proposition 4.3. For any x > and i > 0, it holds 

PALt<i) ~ {S{x)+i)u{{t,+oo)). (4.7) 

Proof. Since Lt increases only when X is at we may write 

P.iLt <i) = P.{Ho >t)+P,{Ho<t,Lt< i) 

= P,{Ho >t) + P,(i7o < t , Lt.Ho o Oho < ^) 
= >t) + P^{Ho<t,t- Ho< n), 

where 9. denotes the usual shift operator and fi is a subordinator starting 
from 0, independent of Hq and identical in law with T£ (under Pq), by the 
strong Markov property. Consequently, 

Px{Lt <i) = P,{Ho >t)+ P,{fe + Ho>t)- P^in + Ho>t,Ho>t) 
= P,if, + Ho > t). 

We use Lemma [231 and take therein F to be the Pa;-distribution fi (which is 
the same as the Po-distribution Tg) and G the Pr-distribution of Hq. Then, 
by (14.21) . F is subexponential and from (14. 3p and ( 14.41) we have 

1„„M^ = ^>0. 

t^^ P,^{Tg >t) I 

Consequently, by Lemma [2?4l 

P,(f, + //o > t) 



lim - 

P^(ff > t) + P^{Hq > t) 



in other words, 



Px{Lt<i) ~ P,{Ho>t)+P,{fe>t) 

~ S{x)u{{t,oo))+iu{{t,oo)), 

t—KX) 

as claimed. □ 

Example 4.4. For a Bessel process of dimension d G (0, 2) reflected at we 
have from (13.281) in Example 13.51 

K(t,+oo)) = ^t-", 
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and Assumption (A) holds by Lemma 12^21 Consequently, 

P,(Lt<£) ~ {S{x)+i)u{it,+oo)). 



t—*oo 



where the scale function is as in Example I3.5[ Taking here a = 1/2 gives 
formulae for reflecting Brownian motion. We remark that our normalization 
of the local time (see (11.21) ) is different from the one used in Roynette et al. 
[22] Section 2. In our case, from f ll.4p and (13.28^ it follows (cf. also [2] p. 
133 where the resolvent kernel is explicitly given) that 

Eo (exp(-Ar,)) = exp (^-i 2^-° A"^ . (4.8) 

Comparing now formula (2.11) in [22] with (14. 8p it is seen that 

Lt = 2aLt 

where L denotes the local time used in [22l. 



5 Penalization of the diffusion with its local 
time 

5.1 General theorem of penalization 

Recall that (C, C, {Ct}) denotes the canonical space of continuous functions, 
and let P be a probability measure defined therein. In the next theorem 
we present the general penalization result which we then specialize to the 
penalization with local time. 

Theorem 5.1. Let {Ft : t > 0} be a stochastic process (so called weight 
process) satisfying 

< E{Ft) < oo V t > 0. 
Suppose that for any u >0 

exists a.s. and 

E(M„) = 1. (5.2) 

Then 
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1) M = {Mu : M > 0} is a non-negative martingale with Mq = 1, 

2) for any m > and A G 

lim = E(UMJ =: Q(")(A), (5.3) 

3) there exits a probability measure Q on {C,C) such that for any u> 

Q(A) = Q(")(A) VA e Cu. 



Proof We have (cf. Roynette et al. [20]) 

E(1a„FJ E(Fi|C.) 



E{Ft) V E(Fi) 

and by (IS.ip and (15.21) the family of random variables 

'E{Ft\Cu 



E(F*: 



: t > 



is uniformly integrable by Sheffe's lemma (see, e.g., Meyer [IT]), and, hence, 
(15.30 holds in L^(f2). To verify the martingale property of M notice that if 
u < V then A^ G Cy and by (15. 3p we have also 

Consequently, 

E(1a„MJ = E(1a„MJ, 

i.e., M is a martingale. Since the family {Q*^"^ : m > 0} of probability 
measures is consistent, claim 3) follows from Kolmogorov's existence theorem 
(see, e.g., Billingsley pj p. 228-230). □ 



5.2 Penalization with local time 

We are interested in analyzing the penalizations of diffusion X with the 
weight process given by 

Ft := h{Lt), t > (5.4) 
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with a suitable function h. In particular, if /i = i[o,e) for some fixed i > then 
Ft = l{Lt<e}- In the next theorem we prove under some assumtions on h the 
validity of the basic penalization hypotheses (|5.1I) and (15 ■2p for the weight 
process {Ft : t > 0}. The explicit form of the corresponding martingale M'^ 
is given. In Section 6.3 it is seen that remains to be a martingale for 
more general functions h, and properties of X under the probability measure 
induced by M'* are discussed. 

In Roynette et al. ^22j this kind of penalizations via local times of Bessel 
processes with dimension parameter d G (0, 2) are studied. Our work gener- 
alizes Theorem 1.5 in |22| for diffusions with subexponential Levy measure. 

Theorem 5.2. Let h : [0, oo) i— > [0, oo) be a Borel measurable, right- continuous 
and non-increasing function with compact support in [0, K] for some given 
K > 0. Assume also that ^ 

/ h{y)dy = l, 
Jo 

and define for a; > 

Then for any u> 

EoihjLt) I 
Boih{Lt)) 

and 



H{x) := / h{y)dy. 



' = S{Xu)h{L^) + 1 - H{L^) =■ a.s. (5.5) 



Eo {M^) = 1. (5.6) 
Consequently, statements 1), 2) and 3) in Theorem \5.1\ hold. 

Proof. J) We prove first (15. 5p . 

a) To begin with, the following result on the behavior of the denominator in 
(15.50 is needed: for any a > 

E,(/i(L,)) ~ (S(a)/i(0) + l) K(t,oo)). (5.7) 

To show this, let /x denote the measure induced by /i, i.e., [Ji{dy) = —dh{y). 
Then 

h{x) = / fx{dy) = / l{y>^^fx{dy), (5.8) 

J(x,K] J{0,K] 
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and, consequently, 

Ea{h{Lt)) = f l{,>Lj/i(rf£)) = [ Pa{Lt < i)fi{de). 

\JiO,K] J J{0,K] 

By Proposition 14.31 

lim —7- = S{a)+ L 

t^oo iy[[t, OO)) 

Moreover, for £ < K 

77 < 77 TT > b[a) + K as t — >• CXD, 

and, by the dominated convergence theorem, 

lim I = / {Sia)+£)M). 



Hence, 



E.(ML,)) ~ f/ iS{a) + i)i^idi)] u{{t,oo)), (5.9) 



'{0,K] 

and the integral in (15. 9p can be evaluated as follows: 



{s{a) + e)fi{de) = s{a) i^{de)+ ifiide) 

{0,K] J{0,K] J{0,K] 

= S{a)h{0)+ [ [ du 

J{0,K] Jo 

= Sia)hiO) + [ du f fx{di) 

Jo J{u,K) 

= S{a)h{0)+ [ h{u)du 
Jo 

= S{a)h{0) + 1. 
This concludes the proof of (15.71) . 

bj To proceed with the proof of (15.51) . recall that {L^ : s > 0} is an additive 
functional, that is, Lt = + Lt-u o 9^ for t > u where 9u is the usual shift 
operator. Hence, by the Markov property, for t > m 

EoihiLt) I Cu) = Eoih{Lu + Lt^u o 9^) \ C„) = H{Xu, L^, t - u) (5.10) 
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with 

H{aJ,r) ■.= Ba{h{i + Lr)). 
By (I5.7p . since x ^— h{i + x) is non-increasing with compact support, 



H{a,i,r)^r^ \^S{a)h{i) + J h{i + u)du j i^{{t, oo)). 

Bringing together ( IS.lOp and (15.7^ with a = yields 

Eo(/i(L,) I C„) _ S{X^)h{Lu) + J^^ h{x)dx 

completing the proof of (15.51) . 

II) To verify (15.61) . we show that {M^ : t > 0} defined in (15. 5p is a non- 
negative martingale with Mq = 1 (cf . Theorem 15.11 statement 1 )) . 
a) First, consider the process S{X) = {S{Xt) : t > 0}. Since the scale 
function is increasing S{X) is a non-negative linear diffusion. Moreover, 
e.g., from Meleard [H], it is, in fact, a sub-martingale with the Doob-Meyer 
decomposition 

SiXt) = Yt + Lt, (5.11) 

where y is a martingale and L is a non-decreasing adapted process. Because 
L increases only when S{X) is at or, equivalently, X is at L is a local 
time of X at 0. Consequently, L is a multiple of L a.s. (for the normalization 
of L, see p. 21) ). i.e., there is a non-random constant c such that for allt > 

Lt = cLt. (5.12) 

We claim that L coincides with L, that is c = 1. To show this, recall that 

Ex(4^^) = / p{s;x,y)ds, 
Jo 

which yields (cf. (II. ip ) 

poo 

Rx{0,0)= / Xe-^'EoiLt)dt. 
Jo 

From ( 15TT11 and (l5T2l l we have Eo(Li) = ^Eo{S{Xt)), and, hence, 

1 f°° 

i?A(0,0) = -/ Xe-^'Eo{SiXt))dt 







- r S{y)XRx{0,y)m{dy). (5.13) 
c Jo 
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Next recall that the resolvent kernel can be expressed as 

Rxix,y) = w^^ipxix)ipxiy), 0<x<y, (5.14) 

where wx is a constant (Wronskian) and ipx and ipx are the fundamental de- 
creasing and increasing, respectively, solutions of the generalized differential 
equation 

d d 

characterized (probabilistically) by 

Consequently, (I5.13P is equivalent with 

1 

1 /"°° d d 

= -cl ^^y^d^Ts'^'^y^'^^'^y^- 

= -/ dS{y) / m{dz)——<^x{z). 
c Jq Jy dm db 

where for the third equality we have used Fubini's theorem. Next we claim 
that 

— cpxi+(x^) := \im— ipxix) = 0. (5.18) 

To prove this, recall that the Wronskian (a constant) is given for all 2; > 
by 

wx = ^z) 4; Mz) + Mz) ( <^a(^) ) . (5.19) 



dS " ' ' ^ ' ' \ dS 

Notice that both terms on the right hand side are non-negative. Since the 
boundary point +00 is assumed to be natural it holds that lim^^oo = +00 
a.s. and, therefore, (cf. (15.161) ) 

lim ipx{z) = +00. 
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Consequently, letting z +00 in (I5.19p we obtain f lS.lSp . Now (15 .17^ takes 
the form 



(^a(0) = — (<^a(+oo)-<^a(0)). 



This implies that c = 1 since (^a(+oo) = by the assumption that +cxd is 
natural (cf. (I51B I 

bj To proceed with the proof that is a martingale, consider first the case 
with continuously differentiable h. Then, applying (IS.lip . 

dM^ = h{Lt){dYt + dLt) + S{Xt)h'{U)dU - h{U)dU = HLt)dYt, (5.20) 

where we have used that 



Consequently, M is a continuous local martingale, and it is a continuous 
martingale if for any T > the process {Mj^ '■ < t < T} is uniformly 
integrable (u.i.). To prove this, we use again (IS.lip and write 



Since h is non-increasing and and has a compact support in [0, K] we have 



showing that {h{Lt)Lt + l — H{Lt) : t > 0} is u.i. Moreover, since {h{Lt) : t > 
0} is bounded and {Yt ■ < t < T} is u.i. it follows that {h{Lt)Yt : < t < 
T} is u.i.. Consequently, {M^ : < t < T} is u.i., as claimed, and, hence, 
{Mj} : t > 0} is a true martingale implying (15. 6p . By the monotone class 
theorem (see, e.g., Meyer |17| T20 p. 28) we can deduce that {M/* : t > 0} 
remains a martingale if the assumtion "h is continuously differentiable" is 
relaxed to be "h is bounded and Borel-measurable". The proof of Theorem 
15.21 is now complete. □ 

Example 5.3. Let h{x) := l[o,f)(x) with £ > 0. Then 



S{Xt)h\Lt)dLt = S{0)h'{Lt)dLt = 0. 



= h{Lt)Yt + h{U)Lt + 1 - H{Lt). 



(5.21) 



\h{Lt)Lt + l 
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and the martingale takes the form 

Mt = j{SiX^)l^L.^e} + ii-Lu)+) 
= J (5(XJ + £ - L„) 1|^„<,} 
~ £ ('S'(Xu/\t-^) + £ — LuAn) 

= 1 + ('S'(X,jAt-^) — LuAtJ ■ 

5.3 The law of X under the penahzed measure 

In this section we study the process X under the penahzed measure Q in- 
troduced in Theorem I5.2[ In fact, we consider a more general situation, and 
assume that h is "only" a Borel measurable and non-negative function defined 
on R_i_ such that 



h{x)dx = 1. (5.22) 







For such a function h we define 

rh 



■= S{Xt)h{U) + 1 - H{Lt), (5.23) 
where 

H{x) ■= [ h{y)dy. 







It can be proved (see Roynette et al. [201 Section 3.2 and [22] Section 3) 
that {Mj} : t > 0} is also in this more general case a martingale such that 
Eo(M/^) = 1 and lim^^oo = 0. Therefore, we may define, for each m > 0, 
a probability measure Q'^ on (C,Cu) by setting 

Q''(A„) := Eo (1a„ M„^) A„ G C„. (5.24) 

The notation E'^ is used for the expectation with respect to Q'^. Next two 
propositions constitute a generalization of Theorem 1.5 in [22] . 

Proposition 5.4. Under Q^, the random variable Loo := linit^oo Lt is finite 
a.s. and 

Q''{L^ G di) = h{i) di. 
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Proof. For m > and ^ > it holds {Lu > i} & Cu, and, consequently, 
By optional stopping, 

Eo (l{Tf<«} = Eo (l{rf<«} AfrJ ) 

but 

M^^ = ^(X.J/i(L.J + 1 - if (L.J = 5(0)/i(£) + 1 - H{£) 

POO 

= / h{y)dy. (5.25) 

As a result, 

Q'^(i:« - = (/ ^^^^ "^^j ^"^""^ - 

Letting here u oo and using the fact that is finite Pq— a.s. shows that 



PC 

Je 



h{y) dy. 



Moreover, from assumption (I5.22P it now follows that L^o is Q^-a.s. finite, 
and the proof is complete. □ 

In the proof of the next proposition we use the process = {Xj : t > 0} 
which is obtained from X (cf. (13.11) ) by conditioning X not to hit 0. The 
process X^ can be described as Doob's /i-transform of X, see, e.g., Salminen, 
Vallois and Yor p3] p. 105. The probability measure and the expectation 
operator associated with X^ are denoted by and E^^, respectively. The 
transition density and the speed measure associated with X^ are given by 

P\t;x,y) := |^^^, m\dy) := S{yf m{dy). (5.26) 



Notice (cf. ([33D) that 



p\t-Q,y) ■.= \imp\t-x,y) = (5.27) 



Consequently, we have the formula 



l = Pj(x/>0)=^ p\t-0,y)m\dy) = fyo{t) S{y) m{dy). (5.28) 
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Proposition 5.5. Let A denote the last exit time from 0, i.e., 

A := sup{t : = 0} 

with A = {•} = 0. Then 

1) Q''(0 < A < oo) = 1, 

2) under Q'* 

a) {Xt ■ t <X\ and {Xxj^t : ^ > 0} are independent, 

b) conditionally on — £, the process {Xt : t < X} is distributed 
as {Xt : t < Tc} under Pq, in other words, 

E"(F(X, : t<\)f{L^)) 

POD 

= / f{i)h{e)Eo {F{Xt -.tKn) dt (5.29) 
io 

where F is a hounded and measurable functional defined in the 
canonical space (C, C, (Cj)) and f : [0,oo) i— [0, cxd) is a bounded 
and measurable function. 

c) the process {Xx+t : t > 0} is distributed as {Xj : t > 0} started 
from 0. 

Proof. Consider for a given T > 

A E'^ (Fi(X„ : K < A) F2{Xx+, : v <T) f{Lx) -'-{0<A<oo}j , 

where F^ and F2 are bounded and measurable functionals defined in the 
canonical space (C, C, (Ct)) and / : [0, 00) 1— > [0, 00) is a bounded and mea- 
surable function. For N > define 

Xn := sup{u <N:Xu = 0} 

and 

aS^ := E'^ ■.u<Xm) F2{Xx,+v : v <T) /(L^,) 1{x^+t<n}) • 

Then 

A = lim aS^\ 
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By absolute continuity, cf. (15.241) . 

Aji^ = Eo (Fi(X„ : u<Xn) F2(Xa^+. : v<T) f{LxJ 
= Eo(Fi(X„ : m < Ajv) F2(Xa,+. ■■v<T) f{Lx^) 

^{Xn+T<N} 

x{SiXN)h{L^) + l-H{L^))y 

Since Fi,F2, and / are bounded and 

lim (1 - H{Ln)) = Po-a.s. 

N—*oo 

we have 

A = lim Eo(Fi(X„ ■.u<Xm) F^iX^^^, : v <T) f{Lx^) 

>^'^{Xm+T<N} S{XN)h{LN)j . 

(2) 

Let A]y denote the expression after the limit sign. Then we write 

A;J) = Eo(^Fi(X„ : u < r,_)F2(X,,_+„ : v<T) 
I 

S{X^)h{l)\, 

where {r^} is the right continuous inverse of {L^} (see (II. 3p ). By the Master 
formula (see Revuz and Yor [18] p. 475 and 483) 

rf£/i(£)/(£)Eo(Fi(X„ : M<r,) 

X y n((ie) F2(e^ : < T) l{T<Af-Tf<c(e)} 'S'(eAr_^J j , 

where E denotes the excursion space, e is a generic excursion, ^(e) is the life 
time of the excursion e, and n is the Ito measure in the excursion space (see, 
e.g., US] p. 480 and p3]). We claim that 

/ := y F^ie-o : v< T) l{T<T'<c(e)} S{eT') n{de) 

= El{F2{X, ■.v<T)). (5.30) 
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Notice that the right hand side of fl5.30p does not depend on T' . We prove 
f Onll for F2 of the form 

F2{ev ■■ V <T) = G{et^,. . .,etj, h < t2 < ■ ■ ■ < tk = T, 

where G is a bounded and measurable function. For simplicity, take k = 2 
and use Theorem 2 in [23] to obtain (for notation and results needed, see 
O, O, iKm and (15:^71) ) 

/xi,o(^i) Pih - h; xi, X2) p{T' - t2] X2, X3) 

[0,oo)3 

X X2) S^Xs) m{dxi) m{dx2) m^dx^) 



S{Xi) fxi,o{tl)pHt2 - ti;Xi,X2)p^{T' -t2;X2,X3) 

[0,oo)3 

X G{xi, X2) S{x2Y S^XsY m{dxi) m{dx2) mi^dx^) 
proving (I5.30p . Consequently, we have (for all N) 

/•OO 

Ag) = El {F2{X, : V < T)) / d^h{^)f{^E, (Fi(X, : u < n)) , 

Jo 

and choosing here Fi, F2, and / appropriately implies all the claims of Propo- 
sition. In particular, Fi = F2 = 1 and / = 1 yields Qo(0 < A < 00) = 1, 
and, hence, L^o = Lx Qg-a.s. □ 



6 Appendix: a technical lemma 

The following lemma could be viewed as a "weak" form of the Tauberian 
theorem (cf. Feller |6] Theorem 1 p. 443) stating, roughly speaking, that if 
two functions behave similarly at zero then their Laplace transforms behave 
similarly at infinity. 

Lemma 6.1. Let be a a-finite measure on [0, +00) and gi and g2 two real 
valued functions such that for some Aq > 

Ci:= [ e-^°^ \gi{^)\ fiid^) < 00, i = l,2. 

J[0,+oo) 



31 



Assume also that (72(7) > for all 7. Introduce for A > Aq 
/,(A):= / e-"^(7.(7)Mc?7), ^ = 1,2. 

"'[0,+oo) 

and suppose 



Then 



implies 



lim /2(A) e"-" = +00 for all 6 > 0. (6.1) 

A— >+oo 



^1(7) ^,^2(7) (6.2) 



/i(A) ^ ~ /2(A) (6.3) 

A— ►+00 



Proof. By property (16.20 there exist two functions 9^ and 9* such that for 
some e > and for all 7 G (0, e) 

9.i6)g2i^)<gi{-f)<9*i6)g2{i). (6.4) 

and 

lime,(£) = limr(£) = 1. (6.5) 

£— >0 e— >0 

We assume also that 9^{e) > and 9*{e) > 0. Letting A > 2Ao we have for 

7 > £ 

A ^ A -^7 ^ A Ae 
A7 > Ao7 + — > Ao7 + Y 

and 



00 



e-"^ Ml)\Kdj) < J g-Ao7 |^^(^)| ;,(ci^) < e-^-/2 Q. (6.6) 

Furthermore, from ( 16. 4p 

e-'^g,{^)f,{d^)<9*{e) f e"^^ (72(7) Mrf7) 

^0 

6-^^^72(7) /^(c^7) 
<^*(^)/2(A) (6.7) 
since g2 is assumed to be positive. Writing 

e-^^^?i(7)/i(rf7) + y e~^^^?i(7)/i(t/7) 
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the estimates in (I6.6p and (16.71) yield 

/i(A)<r(£)/2(A)+e-^^/2Ci, 
which after dividing with /2(A) > implies using (16.ip and (16.5^ 

limsup^^^TT = 1- (6-8) 

For a lower bound, consider 

/i(A)= / e-^^^7i(7)Mrf7) 

^[0,00) 

> / e-^^<7i(7)/i(rf7)- / e-^^|^?i(7)|/i(t/7) 

J[0,e) Je 

>e.{e) [ e-'"' g^ii) fi{d^) - e~'^/' C, 

J [0,6) 

> e.{e) (/2(0) -^"e-^^^?2(7)M^7)) -e'^'^C,. 



Hence, 



/2(A) - ' ' ' '^e^^/2/2(A) 

showing that 

limmf — — > 1, 
/2(A) 

and completing the proof. □ 
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